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ABSTRACT 

We study the correlation between the locations of galaxy-galaxy strong lensing 
candidates and tracers of large-scale structure from both weak lensing or X-ray 
emission. The COSMOS survey is a unique data set, combining deep, high reso- 
lution and contiguous imaging in which strong lenses have been discovered, plus 
unparalleled multiwavelength coverage. To help interpret the COSMOS data, 
we have also produced mock COSMOS strong and weak lensing observations, 
based on ray-tracing through the Millenium simulation. In agreement with the 



1 Laboratoire d'Astrophysique, Ecole Polytechnique Federale de Lausanne (EPFL), Observatoire de 
Sauverny, 1290 Versoix, Switzerland 

2 Astronomisches Rechen-Institut, Zentrum fur Astronomic dcr Univcrsitat Heidelberg, Monchhofstr. 12- 
14, 69120 Heidelberg, Germany 

3 Laboratoire d'Astrophysique de Marseille, CNRS Universite de Provence, 38 rue F. Joliot-Curie, 13388 
Marseille Cedex 13, France 

4 Argclandcr-Institut fur Astronomic, Auf dem Hiigcl 71, 53121 Bonn, Germany 

5 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Strasse 1, D-85741, Garching, Germany 

institute for Astronomy, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, U.K. 

7 Universita "Federico II", Department of Physical Sciences, via Cinthia 9, 80126, Napoli, Italy 

8 Istituto Nazionale Fisica Nucleare - Sezione di Napoli, Italy 

9 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrasse, 85748, Garching, Germany 
10 University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 
n LBNL & BCCP, University of California, Berkeley, CA 94720 

12 Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 

13 Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/SEDI-SAP, Service 
d'Astrophysique, CEA Saclay, Orme des Merisiers, 91191 Gif-sur-Yvette, France 

"California Institute of Technology, MS 105-24, Pasadena, CA,91125, USA 

15 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218 



-2 - 



simulations, we find that strongly lensed images with the largest angular sepa- 
rations are found in the densest regions of the COSMOS field. This is explained 
by a prevalence among the lens population in dense environments of elliptical 
galaxies with high total-to-stellar mass ratios, which can deflect light through 
larger angles. However, we also find that the overall fraction of elliptical galaxies 
with strong gravitational lensing is independent of the local mass density; this 
observation is not true of the simulations, which predict an increasing fraction of 
strong lenses in dense environments. The discrepancy may be a real effect, but 
could also be explained by various limitations of our analysis. For example, our 
visual search of strong lens systems could be incomplete and suffer from selection 
bias; the luminosity function of elliptical galaxies may differ between our real 
and simulated data; or the simplifying assumptions and approximations used in 
our lensing simulations may be inadequate. Work is therefore ongoing. Auto- 
mated searches for strong lens systems will be particularly important in better 
constraining the selection function. 

Subject headings: cosmology, gravitational lensing 

1. Introduction 

Strong galaxy lensing is a powerful tool to probe cosmological parameters. Through 
the angular diameter distances involved in the lensing equation, it is possible to derive the 
Hubble constant (e.g Oguri 2007; Read et al. 2007; Saha et al. 2006) and to put constraints 
on the cosmological constant (e.g. Grillo et al. 2008; Chae 2003). In practice, however, 
such constraints are limited by the accuracy to which the mass distribution of the lens can 
be modeled. One well-known difficulty is the mass sheet degeneracy (e.g. Gorenstein et 
al. 1988; Saha 2000). Moreover, the influence of the environment on the strong lensing effi- 
ciency, and the fraction of lensing galaxies embedded in large scale structure is subject to a 
lively debate (e.g. Treu et al. 2008; Auger 2008; Moller et al. 2007; Meneghetti et al. 2007; 
Oguri et al. 2005; Keeton & Zabludoff 2004; Rusin & Tegmark 2001; Keeton et al. 1997; 
Kochanek 1996). Mass in the environment of a lens galaxy (nearby groups or clusters, or 
intervening mass along the line of sight) is thought to (i) enhance the cross-section of strong 
lensing (Fassnacht et al. 2006; Wambsganss et al. 2005; Keeton & Zabludoff 2004), and (ii) 
increase the angular separation between multiple images (Oguri et al. 2005). Therefore, fail- 
ing to take into account the environment of a modeled lens is likely to bias any cosmological 
results. 

Until today, the main advances in understanding if and how environment influences 
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the formation of strong lenses have come from numerical simulations (e.g. Wambsganss 
et al. 2005; Li et al. 2005, 2006; Hilbert et al. 2007, 2008) and from statistical studies of 
the neighborhood of observed lensed quasars (Moustaka et al. 2007; Williams et al. 2006; 
Oguri 2006; Oguri et al. 2005). Unfortunately, numerical simulations suffer from model as- 
sumptions and finite spatial and mass resolutions; while observations are limited by the small 
number of strong lenses with detailed information about their environment. No convincing 
answer has yet emerged on the qualitative influence of the environment of the galaxy lens 
cross-section. 

The COSMOS project (Scoville et al. 2007) offers a unique opportunity to compare the 
distribution of strong lenses with their environment. This pan-chromatic survey covering 
nearly 2 square degree on the sky has collected the richest data-set in terms of deep multi- 
wavelength observations, providing unique information on the mass distribution in this field. 
In particular, thanks to the COSMOS HST/ACS imaging (Koekemoer et al. 2007) which 
has enabled the measurement of the weak lensing signal of distant galaxies (Leauthaud et 
al. 2007) and photometric redshifts from multi-wavelength ground-based follow-up (Capak 
et al. 2007), Massey et al. (2007) have been able to reconstruct the 3D distribution of mass to 
redshift z ~ 1. Moreover, from the XMM- Newton wide field dataset (Hasinger et al. 2007), 
Finoguenov et al. (2007 and 2008 in prep) have detected about 200 X-ray galaxy clusters 
and groups out to z ~ 1. We have recently (Faure et al. 2008) identified 67 strong lens 
candidates within the COSMOS ACS images: one system of giant arcs in a galaxy cluster 
and 66 strong galaxy-galaxy lenses. The latter include 19 galaxy-galaxy lenses with multiple 
images or strongly bent arcs, unambiguously advocating their lens origin. 

For the first time, we analyze the correlation between the locations of strong lenses and 
large scale structure in the COSMOS field. The results are interpreted in conjunction with 
a simulated survey, obtained from the Millennium Simulation (MS; Springel et al. 2005). 

The paper is organized in the following way. In § [2j we present the strong lens candidates 
sample, the weak-lensing mass-map and the catalog of X-ray clusters and groups used in this 
study. In § [3j we introduce the Millennium Simulation, used in this paper to compare and 
understand the observed distribution of strong lenses in the COSMOS field. The correlation 
between the strong lens sample and the large scale environment is presented in § |4} The 
influence of the environment on the images' angular separation is studied in § |5j In § |6j we 
investigate the correlation between the position of the strong lenses and the X-ray groups 
and clusters detected in the field. The results are discussed in § [7j Throughout this study 
we assume a cosmological mean matter density Qm = 0.25 (in units of the critical density), 
a cosmological constant Q\ = 0.75, a Hubble constant hi o = 0.73 (km s -1 Mpc -1 ), and 
a scale-invariant initial density power spectrum (spectral index n = 1) with normalization 
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cxg = 0.9. Magnitudes are in the AB system. 

2. Data 

2.1. Strong lens sample 

In order to identify strong galaxy-galaxy lenses, Faure et al. (2008) first selected a sample 
9452 of galaxies in the COSMOS-ACS survey, which will henceforth be called the "parent 
population" . This parent population consists of all galaxies in the survey with photometric 
redshifts 0.2 < z p hot < 1.0, absolute V-band magnitude My < —20, and which are classified 
as early-type galaxies according to their color by the photometric-redshift software (Mobasher 
et al. 2007). 

A 10" x 10" region around each galaxy in the parent population was visually inspected 
for multiple images and arc-like features. We also fitted and subtracted the stellar light of 
the lensing galaxy, and built ground-based color images to aid the detection of genuine strong 
lenses. This procedure yielded to 66 galaxy-galaxy lens candidates, which will henceforth 
be referred to as simply "lenses". Among these 66 galaxies, 19 systems are very convincing, 
showing multiple images of similar color or elongated curved arcs. They are called the "best 
systems" throughout this paper. The other 47 systems have generally a single arc, bluer 
than the lensing galaxy, and are referred to as "single arc systems" . 

2.2. Weak lensing mass map 

The weak- lensing (WL) convergence map displayed in Fig. [I] (top) is a measurement 
of the projected mass distribution. The convergence k has been calculated from the image 
shapes of about half a million galaxies detected in the HST/ACS observation (Leauthaud 
et al. 2007, Massey et al. 2007), which have a density ~ 70arcmin -2 and a median redshift 
z~ 1.17. The convergence map has been smoothed using a multi-scale filtering method 
(Starck et al. 2006) with a maximum resolution of 1.2' full- width at half- maximum (FWHM). 
The measurement is sensitive to mass in a broad range of redshift 0.2 < z < 1.0, with a peak 
in the effective sensitivity at z ~ 0.7 (Massey et al. 2007). 

Strongly lensed galaxies can typically lie at slightly higher reshifts (z ~ 2) due to their 
additional magnification. However, the breadth of the lensing sensitivity kernel ensures that 
the weak-lensing mass reconstruction spans almost the same redshift range as the strong lens 
sample. 
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Fig. 1. — Strong lenses projected in the COSMOS weak lensing mass map (top) and in the 
X-ray galaxy clusters and group mass map (bottom). Top: The linear grey-scale indicates 
the convergence (Massey et al. 2007), which is proportional to the (lens-efficiency weighted) 
projected mass along the line of sight. Bottom: the clusters/groups are represented by 2D 
gaussians centered at the location of the cluster/group (FWHM=2xr 50 o, the amplitudes are 
scaled on their weighted surface mass-density; see § [2]). Contours are drawn at k — 0.4% 
and 1.4% in the WL mass-map and reported for comparison in the bottom map. The red 
squares correspond to the single arc systems and the blue squares correspond to the best 
systems. Small squares denote arc radii r arc < 1.5", large squares r arc > 1.5". 
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2.3. X-ray groups and clusters 

Finoguenov et al. (2007, 2008 in prep.) have discovered 202 galaxy groups and clusters 
in the COSMOS field, via X-ray emission from hot gas. This techniques probes virialised 
structures with a mass detection threshold lower than that of the WL mass map, but with a 
different sensitivity as a function of redshift. The galaxy clusters are characterized by a size 
r 500 ( typically ~0. 01°), the radius inside which the matter density is 500 times the critical 
density, and where most of the cluster's X-ray flux is generated (Markevitch 1998). Fig. [I] 
(bottom) displays the X-ray galaxy clusters and groups detected in the COSMOS field. Each 
X-ray source is represented by a 2D gaussian, with FWHM equal to twice the radius r 500 , 
and amplitude scaled to the lensing efficiency of the cluster, assuming background sources at 
z=2. The lensing efficiency is defined by the weighted surface mass-density, £ cx Ms f x ^f^, 

7rXr 500 U S 

where M 50 o is the mass encompassed in the radius r 50 o- 

A comparison of Fig. [ljtop and bottom reveals some discrepancies in the mass distribu- 
tion as traced by weak lensing and X-ray emission. As noted by Massey et al. (2007), this is 
mainly due to the additional sensitivity of the X-ray analysis to compact objects. However, 
the weak lensing map is also noisy and potentially subject to localized defects that are not 
revealed by systematic tests statistically averaged over the entire field. 



3. Simulated data from the Millennium Simulation 

3.1. Ray-tracing setup 

To test whether our observational results are compatible with theoretical predictions, 
we have also created mock COSMOS data by ray-tracing through the Millennium Simulation 
(MS; Springel et al. 2005). This is an A^-body simulation of cosmic structure formation, which 
employed a TreePM version of Gadget-2 with 10 10 particles of mass m p = 8.6 x lO 8 /z _1 M 
to follow the structure formation in a cubic region of L = 500/i -1 Mpc comoving on a side 
from redshift z = 127 and z — 0. 

The ray-tracing through the MS is based on a multiple-lens-plane algorithm described 
in Hilbert et al. (2007, 2008a, 2008b). The dark-matter distribution in the observer's back- 
ward light-cone is generated directly from the particle data of the of the MS, and projected 
onto a series of lens planes. The stellar matter in galaxies is inferred from semi-analytic 
galaxy-formation models implemented within the evolving dark-matter distribution of the 
MS (De Lucia & Blaizot 2007). Light rays are traced back from an observer to their source, 
assuming that the rays propagate unperturbed between lens planes, but are deflected when 
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passing through a plane. The image distortions and amplifications resulting from differen- 
tial deflections of the light rays at the lens planes are calculated from the projected matter 
distribution on the planes. 

3.2. Simulated WL mass maps 

We randomly choose 20 fields of 0.5 x 0.5 deg 2 . In each field, a grid of 8192 x 8192 rays is 
traced through the MS by the multiple-lens-plane algorithm. We calculate the convergence k 
in each field from the ray distortions to sources at redshift z — 1.17 (i.e. the median redshift 
of the galaxies used for the weak-lensing map of the COSMOS field). We then add Gaussian 
noise to each pixel of the resulting convergence maps, with a variance appropriate for a source 
galaxy density of 70arcmin~ 2 and galaxy ellipticity variance ~ 0.3. The noisy convergence 
maps are then smoothed using the multi-scale, wavelet filtering method (Starck et al. 2006) 
that was applied to the real COSMOS map. This COSMOS-like WL simulation will then be 
processed through the same analysis as the real data, to investigate the correlation between 
the strong lenses and large scale structures. 

In Fig. |2j we display the projected mass map in 9 of the 20 simulated fields. They 
contain noticeably less structure than the real COSMOS data. It is already known from a 
study of the angular correlation of the galaxies in the COSMOS field (McCracken et al. 2007) 
that this field is at the upper end of the expected cosmic variance scatter. Hence we expect 
from the numerical simulations to give us qualitative results on the behavior of strong lenses 
rather than direct quantitative results that can be compared to the COSMOS field. 

3.3. Strong-lensing simulations 

From the semi-analytic galaxy catalog of De Lucia & Blaizot (2007), we select all galax- 
ies with redshift 0.2 < z < 1.0, absolute V-band magnitude My < —20, and bulge-to-total 
B-band luminosity ratio -^B,bui g e/-^B,totai > 0.4 (Croton et al. 2006), in order to select similar 
early-type galaxies to those visual inspected in the real data. We calculate the image posi- 
tions of this "parent population" of simulated galaxies by linear interpolation between ray 
positions. 

The identification of strong galaxy-galaxy lenses in the COSMOS field involves the 
human eye, whose "selection function" is difficult to model. In order to model the mag- 
nification bias, we assume an apparent magnitude threshold for detecting strongly lensed 
galaxy images, as well as a redshift and luminosity distribution for sources in the respective 
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filter band. From the faintest arc detected in the Faure et al. (2008) sample, we deduce 
an apparent magnitude threshold of mu m ~ 26 [mag] in the F814W filter band (i.e. close to 
the COSMOS survey F814W magnitude limit; Leauthaud et al. 2007). We assume that the 
number density n s (z s ,mi im ) of sources at redshift z s brighter than the apparent magnitude 
m m i n is simply given by a product 

n a (z s , miim) = Ps(z s , m hm )n s (m hm ) (1) 

of a source redshift distribution 

(3zl exp 



V s {z s ,m llm ) = L 3 — 2 

with parameters /3 = 3/2 and zo( m iim) = 0.13raii m — 2.2, and a cumulative source magnitude 
distribution 

mum 

. (3) 

y'^Q2a(mi-m) _|_ ^Q26(mi-m) 

— oo 

with parameters a = 0.30, b = 0.56, mi = 20, and n = 3 x 10 3 deg~ 2 (Smail et al. 1994, 
Casertano et al. 2000, Leautheaud et al. 2007, Fedeli et al. 2008). Q 

Since most strong galaxy-galaxy lens candidates were selected due to arc-like features, 
we consider the optical depth r for images with length-to- width ratio r > 7 for small circular 
sources (see Hilbert et al. 2008a for details). Furthermore, we only consider strong-lensing 
regions in annuli between 0.2" and 5" around a parent-sample galaxy, since closer images 
could not be distinguished from the center of the lens galaxy, and more distant images are 
outside the visually inspected regions (Faure et al. 2008). 



4. Strong lenses as a function of weak lensing environment 

In this section, we first study whether strong galaxy lenses lie preferentially in dense 
environments, as characterized by a high weak-lensing convergence. Then, by comparison 
with the simulations, we investigate which component of the large structure enhances the 
strong lensing cross-section. 



■"■Note that we also considered a source population with an integrated luminosity function inversely pro- 
portional to the threshold luminosity, and a source population whose images are all detected regardless of 
their brightness. The optical depths for such sources may differ by an order of magnitude, and the ratio 
between the optical depths in different weak-lensing convergence regions may be somewhat lower or higher. 
The qualitative behavior is, however, very similar to the case discussed here. 
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4.1. Correlations in the real data 

We split the WL mass map of the COSMOS field in three regions according to their 
convergence values: a low-density region with k < 0.4%, an intermediate-density region with 
0.4% < k < 1.4%, and the high-density region with k > 1.4% (see Fig. [I]). The limits 
between the regions are set arbitrarily, and are only aimed at studying a possible transition 
in the strong lenses behavior of galaxies between regions of low and high convergence. The 
maximum value measured in the WL map is k ~4.5%. Such a low maximum convergence is 
due to the finite spatial resolution (>2') of the WL mass map reconstruction, which washes 
out the stronger signal from the core of centrally concentrated matter structures. For the 
same reason, the convergence measured in the WL mass map cannot be accounted for the 
convergence that directly affects the strong lenses, as what is of matter for strong lensing is 
the fluctuation that is coherent over the angular scale of the lens (i.e ~1"). For each of the 
three defined convergence regions, Table [T] lists the density of elliptical galaxies and strong 
lenses, as well as the fraction of galaxies from the parent population that are strong lenses, 
best systems or single arc systems. The error bars correspond to simple Poisson statistics. 

The results are summarized in Fig. [3j We find a comparable density of strong lenses in 
the low mass density environment (all lenses: 35 ± 6 deg -2 , best systems: 12 ± 3 deg~ 2 ) and 
in the high mass density environment (all lenses: 45 ±22 deg -2 , best systems: 22 ±16 deg -2 ). 
If we take into account the fact that the density of elliptical galaxies is about twice as high in 
the densest regions of the field than in the low-density region, we also find a similar fraction 
of elliptical galaxies that are lensing galaxies in high-density regions (all lenses: 5.4 ± 2.7%o, 
best systems: 2.7 ± 1.9%o) and in the low density regions (all lenses: 7.9 ± 1.2%o, best 
systems: 2.6±0.7%o). 

Surprisingly, we thus do not observe any particular correlation between the environment 
and the frequency of strong lenses. Strong lenses appear to be equally distributed throughout 
the different weak-lensing convergence regions of the field. 

4.2. Correlations in the simulated data 

We now analyse our simulated data, to understand whether the absence of observed 
correlation between strong galaxy-galaxy lensing and weak-lensing convergence is compat- 
ible with the standard model of structure formation and galaxy evolution. The lensing 
simulations may also help to interpret the COSMOS observations. 

The results of our lensing simulation are summarized in Table [2] We use the same 
threshold values for the convergence (i.e. k =0.004 and 0.014) in the smoothed convergence 
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Table 1: Summary of the measurements in the weak-lensing mass map. Columns 1-3: The 
three convergence regions described in § |4.1[ Line 1: Covered area. Line 2: Covered area 



fraction. Lines 3-8: Number of galaxies in the parent catalog, number of all strong lens 
candidates, single arc systems, best systems, r^™ a11 systems and r^ ge systems. Lines 9-14: 
Density of galaxies, strong lenses, single arc systems, best systems, and r a ™ a11 and r^ ge 
systems. Lines 15-17: Fraction of galaxies from the parent catalog which are strong lenses, 
single arc and best systems. 
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r arc > 1-5" systems 
Galaxy density [deg -2 ] 

Parent catalog 

All lenses 

Single arc systems 

Best systems 
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Table 2: Summary of the results from ray-tracing of 20 fields of 0.5 x 0.5 deg 2 through the 
MS. Given are the values obtained from the area of all fields combined with errors indicating 
the maximal deviation observed in individual fields. Line 1: Area fraction. Line 2: Number 
density n of galaxies in the simulated parent catalog. Lines 3-5: Strong- lensing optical depths 
r for various source redshifts z s . Lines 6-11: As above, but restricted to systems with arc 
radius r arc < 1.5" and r arc > 1.5", respectively. Lines 12-14: Strong-lensing cross-sections 
a for various source redshifts z s . Line 15: expected number N L of strong lenses (without 
Poisson errors) assuming a total survey area of 1.8 deg 2 . Line 16: expected fraction of strong 
lenses among the parent population (without Poisson errors). 
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maps to divide our simulated lensing fields into low-, intermediate-, and high-convergence 
regions as for the observed field. The resulting area fraction covered by the low- (75 ± 9%), 
intermediate- (19 ± 7%), and high-convergence regions (7 ± 3%) in the simulated fields is 
very similar to the fractions measured in the COSMOS field (except that the fraction of the 
simulted survey area with high weak-lensing convergence is somewhat smaller than that in 



The galaxy densities n of the simulated parent population in the low- and intermediate- 
convergence regions are slightly smaller than the densities observed in the COSMOS field. 
However, except for the low-convergence region, the observed densities are within the broad 
range spanned by the values calculated in individual ray-tracing fields. We are therefore 
confident that our lensing simulations reproduce the properties of the lensing observations 
reasonably well for a quantitative comparison of the correlation between strong galaxy-galaxy 
lensing and weak-lensing convergence. 

As a simple approximation, we assume that each lens galaxy induces one strongly dis- 
torted image (a "single arc") of a source galaxy. Thus, the density n L of strong-lens galaxies 
is given by an integral 



of the optical depth r(z s ,mi im ) and the source density n s (2 s ,mii m ) for given image-detection 
limit magnitude mr m over all source redshifts z s . In the following, we assume m ^ m = 26. 

We are primarily interested in the relative number of expected strong galaxy-galaxy 
lenses in the different convergence regimes. Therefore, we first discuss the optical depth and 
take the source redshift distribution and density into account later. The optical depth t(z s ) 
in the low, intermediate and high weak-lensing convergence regions is shown in Fig. [1] (a) 
for source redshifts z s — 1, z g — 2, and z s = 3. There is a clear trend of an increas- 
ing strong-lensing optical depth with increasing weak-lensing convergence for all considered 
source redshifts. On average, the optical depth in the high-convergence region is 9-15 times 
larger (depending on the source redshift considered) than in the low-convergence region. This 
strong increase is in contrast to the weak increase of the density of strong-lens candidates 
seen in the COSMOS field. 

The increased strong-lensing optical depth r in regions of higher weak-lensing conver- 
gence is partly due to the increased galaxy density n of the parent population. However, 
an increase in the strong-lensing probability remains even if the effect of the larger galaxy 
density is taken out. This can be seen in Fig. [1] (b), where the average strong-lensing cross- 
section a = r/n of galaxies in the parent population is shown. The average cross-section in 
high-convergence regions is 3-7 times larger than in low-convergence regions. Thus, the ex- 



the COSMOS field). 




(4) 
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pected fraction of strong lenses within the parent population in regions with high convergence 
is roughly 3-7 times larger on average than in regions of low convergence. Among individual 
ray-tracing fields, the ratio of cross-sections in the low- and high-convergence region ranges 
from 2 to 10. Hence, the simulations suggest that one should expect the fraction of strong 
lenses among the parent population to be at least 50 percent higher in the high-convergence 
regions than in the low-convergence regions (and even much larger ratios of 3-7 are expected 
on average). 

The expected number of strong lenses in the different weak-lensing convergence regions 
is given by the product = An^ of the area A and lens density in the considered region. 
The expected fraction of galaxies in the parent population that are strong lenses is given 
by the ratio n^/n^. The resulting strong-lens numbers, assuming a total survey area of 
1.8 deg 2 , and the strong-lens fraction are shown in Table [2} For the high-convergence region, 
the predicted number and fraction of strong lenses are consistent with the values observed in 
the COSMOS field. However, the simulation predicts far fewer galaxy-galaxy lenses in the 
intermediate- and low-convergence region than there are lens candidates in the COSMOS 
field. 

The reasons for the different discrepancies between observations and simulations are dis- 



cussed in 5 4.3 



4.3. Observations versus simulations 



From Tables [2] and [TJ we see that the surface covered by the three convergence regions 
as well as the distribution of parent population galaxies is comparable in the simulated 
fields and in the observations. However, the increasing fraction of strong lenses in regions 
of increasing convergence expected from the simulations is not matched by the constant 
fraction of observed strong lenses in the COSMOS field. 

The discrepancy between the simulations and observations may be due to a number of 
reasons. Some of the possible reasons related to the observations are: 

(i) The sample of strong lenses might be incomplete in a way that biases the results. 

(ii) A significant fraction of the strong-lens candidates are not strong galaxy-galaxy lenses. 

(iii) The measured WL mass map does not represent the actual mass distribution in the 
COSMOS field very accurately. 

(iv) Our arbitrary cuts in WL convergence might not encompass important transitions 
between low- and high-density regimes. 
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(v) The COSMOS field might be a peculiar case due to cosmic variance. 

(vi) The number of strong lenses might be too low to reach the necessary statistical preci- 
sion, and therefore we are missing the correlation. 

(vii) There may be actually no correlation between the strong lenses and the environment 
(at least for the used galaxy sample selection criteria). 

Point (i) cannot yet be completely ruled out. Fig. [5] (upper plot) shows the distribution 
of arc radii in our strong lenses. Note the conspicuous absence of observed systems with 
arc radii around ~1" or greater than 3.5". A much smoother distribution is seen in the 
simulations: we would have expected to find systems with arc radii up to 5", which is still 
well within the 10" x 10" visually inspected regions. The gap at 1" corresponds suspiciously 
to the mean seeing size of the ground based imaging. It is possible that genuine strong 
lenses were excluded during the transition from monochromatic, space-based images used to 
detect potential candidates, to ground-based, color images used to verify them. However, 
a second independent search (Jackson 2008) did not discover any new strong lens systems 
amongst the Faure et al. (2008) elliptical galaxies, so we are confident that our sample is 
close to complete. Confirming the statistical significance of the apparent lack of lenses with 
arc radii above 3.5" will require a larger sample. Overall, if the properties of strong lenses 
were strongly influenced by the environment of the main lensing potential, we would have 
expected a much smoother transition between systems with a galaxy-like (r arc <4") and 
cluster-like (r arc >10") mass deflectors, as seen in the simulations (Fig. |5j bottom plot). 

Point (ii) is currently being addressed with follow-up observations (Faure et al. in prepa- 
ration). However, explaining this discrepancy would require a strong bias in the selection 
process as a function of environment - for example favoring the spurious detection of false 
positives in low density environments. Point (iii) is also subject to ongoing methodological 
improvements. The insignificant reconstructed I?-modes (Massey et al. 2007) demonstrate 
the overall reliability of the weak lensing map. However, isolated imperfections in the anal- 
ysis (e.g. imperfect PSF modeling in one or two HST/ACS pointings) can go undetected 
when averaged over the entire field. Such artefacts (as well as noise) can manifest as spu- 
rious convergence peaks, which would dilute the fraction of strong lenses measured in high 
density environments. It is difficult for systematics to remove real lensing convergence and 
thus produce the opposite effect; the much larger area of low density than high density also 
means that noise acts in the same sense. This bias doubtless contributes to the observed 
results, but is not at a sufficient level to explain them alone. Point (iv) can be excluded 
if the lensing simulations are not terribly unrealistic, since they show a clear change in the 
strong-lens fraction. Point (v) can also be excluded, the measured scatter between the 20 
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simulated fields showing that cosmic variance alone is insufficient to explain the discrepancy. 
Future observations of larger fields will also test this. They will also be required to test 
explanations (vi) and (vii), which cannot be ruled out with existing data. 

Apart from possible reasons related to the observations, the discrepancy may be caused 
by the simple assumptions and approximations used in the lensing simulations. These include 
the assumption of circularly symmetric projected stellar-mass profiles of the lens galaxies (i.e. 
lensing cross-section), the neglect of the response of their dark-matter halo to baryon cooling 
(Gnedin et al. 2004), the assumption of point-like sources (instead of extended elliptical 
sources), and a simple source magnitude distribution (which neglects any evolution with 
redshift). In principle, we also cannot exclude the possibility that the discrepancy is due to 
shortcomings of the underlying dark-matter structure formation model. Moreover, the model 
for the strong-lens selection function of the COSMOS sample that we use in the simulations 
might be unrealistic. Firstly, the simulated parent population relies on semi-analytic models 
to predict the galaxy properties; the number of small galaxies in high-density regions could 
have been overestimated. Then the detection criteria for strongly distorted images may be 
too simplistic. For example, the uniform magnitude detection threshold does not depend on 
the distance to the lens galaxy (apart from a cutoff for images closer than 0.2"), or on the lens 
galaxy's apparent size and brightness. In the real data, faint lensed images at small projected 
distances from a large and bright lensing galaxy may have been missed. In addition, there is 
considerable scatter in the average cross-sections measured in different simulated fields (see 
Table [2]). A different value for the normalization erg might change the total number of strong 
lenses, but not their distribution among the different density region. 

In summary, there remain many questions (which we intend to address in future work) 
before we can firmly conclude whether the observed lack of correlation between strong lenses 
and their environment is real, or whether predictions from current models are right. In 
particular, a sample of strong lenses identified via an automated detection algorithm would 
have a better constrained selection function. This would allow further discussion of the 
potential inadequacies of the simulations. 

4.4. The properties of the lens galaxies 

The semi-analytic galaxy catalogue (De Lucia & Blaizot 2007) used for the lensing sim- 
ulations allow us to investigate the relation of various galaxy properties to observed trends 
in the strong lensing probabilities. Fig. |6](a) shows the distribution n~ 1 dn/d\ogM* of the 
stellar mass M* in galaxies of the parent population, for the three convergence regions. 
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The distributions are very similar apart from a small shift towards higher stellar masses for 
regions with larger weak-lensing convergence. This shift, however, implies a higher num- 
ber of galaxies with stellar masses M* > 10 /t M . Although a small fraction of the 
entire parent populations (~4% of the parent population in low- convergence regions, ~9% 
in high-convergence regions) these galaxies contribute significantly to the average lensing 
cross-section of galaxies, as seen in Fig. ^h). These massive galaxies not only have higher 
stellar masses (and thus higher luminosities), but are often associated with significantly more 
massive dark-matter halos associated with them, resulting in very large total-to-stellar-mass 
ratios M tota i/M s teiiar > 100 (and large mass-to-light ratios), as can be seen in Fig. [7j There 
are only a few of them, so they do not significantly affect the luminosity distribution of the 
parent population, but they dominate the lensing statistics. 

This result of the simulation is not in contradiction with the observations in the COS- 
MOS field: the brightest galaxies do not dominate the population of galaxies in the high- 
convergence region (Fig. [8]). To verify the existence of these high total-to-stellar mass galax- 
ies, galaxy-galaxy- weak lensing analysis of the COSMOS galaxies will be necessary (such as 
Mandelbaum et al. 2005), as well as a detail mass model of the strong lenses of the sample. 
Nevertheless, the low correlation between the strong lenses and the mass map measured 
in the COSMOS field suggests that there are fewer of theses massive galaxies in the ob- 
served field rather than in the simulated fields. This result, if verified, could also pinpoint a 
problem in the semi-analytic model used to represent the galaxy luminosity function in the 
simulations. 

The higher cross-section in the high-convergence regions predicted by the simulations is 
not only due to the higher number of massive galaxies. In Fig. [9j the (average) strong-lensing 
cross-section a(M*) = (dr/dM*)/(dn/dM*) as a function of stellar mass M* is shown for 
the different weak-lensing convergence regions. Even if we compare galaxies at fixed stellar 
mass M*, the cross-section a(M*) is consistently larger in the high-convergence regions than 
in the low-convergence regions. Using the simulation data, we can identify at least one 
cause: For given stellar mass, galaxies in high-convergence regions tend to have larger dark- 
matter halos than galaxies in low-convergence regions. However, other possible causes, e.g. 
additional matter from group/cluster halos or other structures along the line-of-sight, cannot 
be excluded so far (and should be investigated in more detail in future work). 

5. Arc radii as a function of weak lensing environment 

The images of a source deflected by a gravitational lens are formed at stationary points of 
the arrival-time surface. The image positions are therefore directly related to the convergence 



-17- 



and the associated shear of the total lens potential. One may conjecture that lens galaxies in 
high-convergence regions will generate larger image separations for multi-image systems than 
lens galaxies in low-convergence regions. Using a model derived from iV-body simulations and 
galaxy formation models, Oguri et al. (2005) studied the impact of external convergence and 
shear on the image-separation distribution of quasars that are multiply imaged by galaxies. 
They found that the angular separation between multiple images is strongly correlated with 
the surface mass density of the environment and that the external convergence enhances the 
lensing probability by ~ 30% for systems with arc radius of r arc « 1.5" (and up to 200% for 
r arc w 2.5"). 



5.1. Arc radii in the COSMOS field 

Here we want to verify whether the distribution of strong galaxy-galaxy lenses arc radii 
in the COSMOS field is a function of the environmental convergence. For that purpose, 
we split the lens sample into two (similarly sized) sub-samples according to their arc radii: 
r|™ a11 if r arc <1.5" and r^ ge if r arc > 1.5". The largest angular separation in our sample is 
r arc = 3.54". The number and density of strong lenses as a function of environment are listed 
for the two sub-samples in Table [TJ and their densities are plotted in Fig. 10 We find that the 



densities of r|™ a11 and r^ ge systems are similar in low- and intermediate-convergence regions. 
However, only r^ ge systems are found in high-convergence environments (48±24deg -2 ) and 
their density is up to four times higher than in low-convergence regions (15 ± 3deg -2 ). 



5.2. Arc radii in the simulated fields 

To estimate the density of strong-lens systems with small and large arc radii, we sepa- 
rately calculate the strong- lensing optical depth for systems with arc radii r arc < 1.5 " and for 



systems with r arc > 1.5 ". The resulting optical depths are shown in Fig. 11 and given in Tab 
[2] for different source redshifts. The optical depth for r^p systems increases very strongly 
with increasing weak-lensing convergence (~20 times higher in high-convergence regions). 
This matches the observed data, and is due to a combination of external shear /convergence 
plus an increased fraction of high mass-to-light ratio galaxies in dense environments. 

The optical depth for r|™ a11 systems also increases (~ 5 times higher in high-convergence 
regions) in the simulations. An equivalent increase in the number of observed r^. systems 



is not seen in the COSMOS data. However, as already mentioned in § 4.3, faint lensed 



images at small projected distances from a large and bright lens galaxy may not be visible in 
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real images, thereby leading to an underestimate of the number of small arc radius systems. 



5.3. Conclusions on the arc radius distribution 

High-density environments enhance the production of strong lens systems with large 
image separations with respect to low-density environments. This effect is observed both 
in the simulations and in the COSMOS field. The large image separations are due to the 
presence of high mass-to-light ratio galaxies in the dense regions and the strong external 
shear /convergence generated by the lens-galaxy environment. The latter can also be induced 
by massive structures along the same line-of- sight. 



As discussed in § 4.4 we do not find substantial evidence from the galaxy luminosity 
function (Fig. [8]) that the brightest galaxies (also likely to be the most massive) are prefer- 
entially found in high density environments. Nor does the low (constant) fraction of strong 
lenses in high-convergence regions of the COSMOS field hint at a excess of massive galaxies. 
However, the presence of only strong lenses with large image separations implies that this is 
indeed the case; or simply that it is only external shear that increases the image separation, 
not the presence of massive galaxies. 



6. Strong lenses as a function of X-ray environment 

As is apparent in Fig. [TJ the WL mass map and the A— ray galaxy clusters and groups 
do not trace exactly the same mass distribution, the latter showing peaks of mass with low 
spatial extension missed by the WL mass map. In order to probe the correlation of strong 
lenses with their environment at higher resolution than the weak lensing mass map, we have 
studied the correlation between strong lenses and their environment as traced by A— ray 
emission (see § 2.3). 



6.1. Lensing galaxies lying in the cluster/group center 

Seven lensing galaxies (2 "best system" and 5 "single arc" systems) are members of 
an A— ray galaxy cluster or group (lying inside a projected radius r 50 o and within redshift 
Az = ±0.05). These represent 10±4% of the strong lenses (10±7% of the "best systems"). 

Two lensing galaxies (1 "best system" and 1 "single arc"; both r^ ge systems) are within 
only 0.02 r 500 of the cluster core. The other five lensing galaxies (of which two are r^ ge 
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systems) lie at projected distances > 0.3 r 500 from the cluster center. 



6.2. Lensing galaxies aligned with a galaxy cluster 

We find that 14 systems (6 "best systems" and 8 "single arc" systems) have an impact 
parameter with a galaxy cluster smaller than rsoo, but have a different redshift than the 
cluster. By adding the 7 systems embedded in a galaxy cluster, we find that 31±7% of 
strong lenses (and 42±15% of the best systems) are within or aligned with a galaxy cluster, 
with an impact parameter inside r 50 o- In comparison, Williams et al. (2006), in a study of 
the environment of 12 strongly lensed quasars with similar arc radii to our sample, found 
that 50±20% to 67±23% of strong lenses are aligned with a galaxy cluster /group. Within 
uncertainties, we are thus in rough agreement with their results. 

If we assume that galaxy clusters in the COSMOS field are circular, do not overlap and 
have a projected size of 7rr| 00 , they cover a projected area of 0.190 deg 2 . We find a density 
of 110 ± 24deg~ 2 strong lenses inside this area, and a density of 28 ± 4 deg -2 outside it 



see Fig. 12). The density of strong lenses is thus ~4 times higher in the projected region 



covered by galaxy clusters. However, the density of galaxies from the parent population 



(Fig. 12) is also ~ 4 times higher in this region. The fraction of elliptical galaxies that are 
strong lenses is therefore approximately constant, regardless of whether the lensing galaxy is 
enclosed within (6±l%o) a projected radius r 50 o from the center of a galaxy cluster or groups 
or not (7±l%o). This recovers the result of §|2J that the population of strong lenses follows 
the distribution of bright elliptical galaxies, and that additional matter structures along the 
line-of-sight do not significantly increase the strong lensing cross-section. 

The density of r|™ a11 and r^ ge systems is similar in the regions away from lines of sight to 
clusters (14 ± 3 deg -2 and 11 ± 3 deg -2 respectively). However, the density of r^ ge systems 
is slightly higher (63 ± 18 deg -2 for r^ ge and 47 ± 16 deg -2 for r^™ a11 systems) when the 
projected distance of a lensing galaxy to a galaxy cluster center is within r 500 . As in § [5j 
two explanations are possible: additional matter along the line of sight increases the angular 
separation between the images, and/or more massive galaxies are found in clusters, which 
boost the lensing cross-section. 



Fig.[T3]shows the distribution of elliptical galaxies in the parent population, according to 
their brightness and as a function of their location. We find a similar distribution of galaxies 
inside and outside the region covered by clusters. An excess of bright, massive galaxies is 
therefore not a plausible explanation for the excessive fraction of r^ ge in the region covered 
by the clusters, unless the galaxies in this region have an intrinsically higher mass-to-light 
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ratio than those outside. 



7. Summary and discussion 

We reach the same conclusion by studying the correlation between strong gravitational 
lenses and large scale structure, as traced by both weak-lensing convergence and X— ray 
emission. We find that dense environments do not enhance the fraction of galaxies that act 
as strong lenses (with image separations r arc <4"), but do increase the number of lenses with 
large image separations. 

We have constructed mock observations of gravitational lensing by ray-tracing through 
the Millennium Simulation (Springel et al. 2007). Simulated regions of high weak-lensing 
convergence similarly increase the number of strong lenses with large image separations. 
However, in contrast to our observations, the simulations also predict a higher total fraction of 
strongly lensing galaxies in high convergence regions than in low convergence regions. Several 



possible reasons for this discrepancy are discussed in § |4.3[ including imperfect modeling 
of the parent galaxy sample, or possible bias in the strong lens selection function. For 
example, the apparent magnitude and redshift distribution of the source galaxy population, 
in particular for high redshifts and faint magnitudes, should be improved, since these are 
crucial for modeling the magnification bias in strong lensing observations. Nor did we take 
into account the shapes and sizes of sources, images or lens galaxies, and may have thus 
missed real strongly lensed images close the lens center (Meneghetti et al. 2008). In future 
surveys, automatic detection software, such as Arcfinder (Seidel & Bartlemann 2007) or 
Haggles (Marshall et al. 2008), will be important to quantify the selection bias. This 
should ease the comparison between observations and simulations. We will investigate such 
possibilities in a further paper (Faure et al. in preparation). 

According to the Millennium Simulations, the main cause for the higher fraction of 
strong lens systems with large image separations in high-density regions is an excess of 
galaxies with high total-to-stellar mass ratio. However, neither the existence of such galaxies, 
nor the fraction of galaxies that they constitute within the field is entirely clear. Detailed 
mass modeling of the various galaxy populations via galaxy-galaxy weak lensing analysis 
(Leauthaud et al. 2008 in preparation) will help explore the existence, location and lensing 
properties of such massive galaxies. 

Constraining simulations and cosmological parameters will ultimately require a larger 
sample of strong lenses. However, to understand the effect of local environment, it is vital 
that the environment be well measured: a task for which the COSMOS survey has been 
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uniquely well designed. The CFHTLS Strong Lensing Legacy Survey (SL2S, Cabanac et 
al. 2007) also aims to explore the relation of strong-lenses to their environment, but is only 
sensitive to larger Einstein radius systems. The final word may only come from the proposed 
SNAP/JDEM mission, which will probe a full range of arc radii. 

We are thankful to the referee for its useful report. JPK acknowledges support from 
CNRS, CNES and the ANR through the grant 06-BLAN-0067. SH is supported by the 
DFG within the Priority Programme 1177 under projects SCHN 342/6 and WH6/3. RM is 
supported by STFC Advanced Fellowship PP/E006450/1. 
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Fig. 2. — Weak-lensing mass maps from the Millennium Simulation. Each map is 0.5° x 0.5°. 
The green squares show the location of the strong lenses. Small squares denote arc radii 
r arc < 1-5", large squares r arc > 1.5". Contours are drawn at k — 0.4% and 1.4%. 
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Fig. 3. — Density of strong lens candidates as a function of the convergence region. The 
density of the parent population has been divided by 100 in this graph to allow comparison. 
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Fig. 4. — Strong-lensing optical depth r and cross-sections a in the low, intermediate and 
high weak-lensing convergence regions for sources at various redshifts z s . The error bars are 
the dispersion of the values measured in individual fields. 
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Fig. 5. — Distribution of arc radii in the COSMOS field (upper plot) and in the simulated 
fields (bottom plot). [Upper plot]: all the COSMOS lenses (solid lines), only the best 
systems (dotted lines). Note the lack of "best systems" at r arc ~ f arcsec. [Bottom plot]: 
for the whole sky (solid lines), for the high-convergence regions (dotted lines) and for the 
low convergence regions (dashed lines). 
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Fig. 6. — Stellar-mass distribution 7i _1 dn/dlogM* of the parent population and the dis- 
tribution of the average lensing cross section of the parent population n _1 dT/dlogM* (for 
sources at redshift z s = 2) as a function of the stellar mass M* for the different weak-lensing 
convergence regions. 
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Fig. 7. — Distribution of the stellar mass M* and total mass M tot of the parent population 
for the low-convergence regions (a) and the high-convergence regions (b). 
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Fig. 8. — Distribution of elliptical galaxies from the parent catalog versus their absolute 
magnitude in the V-band. Top panel: Solid line histogram: galaxies in the low+intermediate 
convergence regions. Dashed line: galaxies in the high convergence level region. Error bars 
are poisonnian. Bottom panel: Ratio between the histograms values. 




Fig. 9. — Cross-section a(M*) of galaxies with stellar mass M* for strong lensing of sources 
at redshift z s = 2 for the different weak-lensing convergence regions. 
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Fig. 10. — Density of strong lens candidates (in deg 2 ) in the three convergence regions 
according to the arc radius. 
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Fig. 11. — Strong-lensing optical depth r for sources at redshift z s = 2 in the low, inter- 
mediate and high weak-lensing convergence regions for systems with arc radius r arc < 1.5" 
and for systems with r arc > 1.5". The error bars indicate the observed range of values in 
individual fields. 
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Fig. 12. — Density of strong lenses off-center with the central region of a cluster (outside 
r 500 ) and aligned with the central region of a galaxy cluster /group, or in it at a distance 
r 500 < 1- The density of the parent population has been divided by 100 in this graph to 
allow comparison. 



Fig. 13. — Distribution of elliptical galaxies from the parent population versus their absolute- 
magnitude in the V-band. Top panel: Dashed line histogram: the distribution of galaxies in 
the region covered by the galaxy clusters (in t^q). Solid line histograms: the distribution of 
galaxies outside this region. Bottom panel: ratio of these two distributions 



